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EQUILIBRIUM  AND  WAVE  PROPERTIES  OF  TWO-DIMENSIONAL  ION  PLASMAS 


Recent  experiments  have  demonstrated  the  existence  of  collective  modes  in 
a  two-dimensional  single  component  ion  layer  residing  just  below  the 
surface  of  liquid  helium.  A  theory  of  the  equilibrium  and  the  wave  properties 
of  such  a  system  is  presented.  The  equilibrium  is  calculated  by  balancing  the 
repulsive  self-electric  field  pressure  with  the  confining  external 
electric  field.  In  the  limit  of  temperature  T  =  0,  the  plasma  is  in  the  form 
of  a  constant  density  disk  at  the  edge  of  which  the  density  decreases  to  zero 
with  a  scale  length  set  by  the  dimensions  of  the  enclosing  cylindrical  box; 
increasing  T  increases  the  width  of  the  transition  region.  Modelling  the 
ions  as  a  cold  two-dimensional  fluid,  it  is  found  that  the  plasma  supports  an 


infinite  set  of  radial  modes  for  each  value  of  (*,  the  azimuthal  mode  number. 
Imposition  of  a  constant  magnetic  field  perpendicular  to  the  charge  sheet 

i  C- 

increases  the  frequency  of  the  1=0  modes  and  for  Q,  splits  each 

mode  into  two;  the  lowest  of  these  split  modes  is  related  to  the  diocotron 


mode 
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I.  INTRODUCTION 

Recent  experiments*  have  succeeded  in  confining  helium  ions  under  the 
surface  of  liquid  helium.  The  ions  form  a  nearly  ideal  two-dimensional  charge 
disk  held  in  position  by  static  electric  potentials  applied  to  the  walls  of 
the  confining  pillbox-shaped  cell  and  by  electric  fields  due  to  the  polariza¬ 
tion  of  liquid  helium  (see  Section  II).  Waves  can  be  excited  in  the  static 
ion  system  by  applying  an  oscillating  potential  to  a  wall  of  the  cell  and 
sharp  resonances  at  discrete  values  of  the  frequency  have  been  observed. 

A  theory  of  the  equilibrium  and  the  linear  wave  properties  of  the  ion 
system  is  presented  here.  The  theory  is  based  on  the  model  of  a  classical 
two-dimensional  fluid  plasma.  The  equilibrium  is  calculated  by  balancing  the 
repulsive  self-electric  field  and  the  plasma  pressure  by  the  confining 
external  electric  field.  In  the  limit  of  zero  temperature  it  is  found  that 
the  self-consistent  plasma  density  is  almost  constant,  out  to  some  radius 
where  the  density  falls  to  zero  on  a  scale  set  primarily  by  the  height 
of  the  confining  cell.  When  the  temperature  is  not  zero,  a  Debye  sheath  is 
added  to  the  transition  region,  making  it  wider  for  increasing  values  of  the 
temperature.  For  typical  experimental  parameters,  the  ratio  of  the  effective 
Debye  length  to  the  height  of  the  cell  is  i  3  *  10-^  and  hence  the  zero  tempera¬ 
ture  density  profile  is  an  excellent  approximation  to  the  actual  profile. 

To  investigate  the  small  amplitude  electrostatic  plasma  oscillations  of 
the  ion  disk,  the  equilibrium  density  profile  in  the  zero  temperature  limit  is 
used  in  the  linearized  fluid  equations  (governing  the  two-dimensional  motion 
of  the  ion  plasma)  and  in  Poisson’s  equation  (which  relates  the  three- 
dimensional  potential  to  the  ion  charge  density).  Normal  modes  of  the  system 
are  determined  by  imposing  the  condition  that  the  wave  potential  vanishes  on 
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the  conducting  wall  of  the  confining  cell.  It  is  found  that  an  infinite  sec 
of  radial  modes  is  associated  with  each  azimuthal  mode  number  l.  The  experi¬ 
mentally  observed  resonances2  can  be  identified  with  the  theoretical  modes 
corresponding  to  low  values  of  the  radial  and  azimuthal  mode  numbers. 

An  interesting  consequence  of  the  polarizability  of  liquid  helium  shows 
up  in  the  plasma  wave  properties.  One  may  at  first  sight  expect  that  the 
plasma  particle  mass  (which  enters  into  the  determination  of  the  resonant 
frequencies)  is  ~  4  a.u.  corresponding  to  an  He+  ion.  However,  the  force  of 
attraction  between  an  He+  ion  and  the  surrounding  polarized  atoms  of  the 
liquid  is  estimated  to  provide  enough  pressure  to  solidify  around  each  ion  a 
casing  of  -  20  atoms  which  moves  together  with  the  ion.  Furthermore, 
since  this  large  particle  moves  in  a  fluid,  its  effective  mass  is  further 
enhanced3  by  half  the  mass  of  the  displaced  fluid.  Thus,  the  effective  mass 
is  expected  to  be  **  120  a.u..  The  value  116  a.  u.  gives  the  best  agreement 
between  theory  and  experiment2.  This  large  value  of  the  effective  mass  has  the 
practical  advantage  of  lowering  the  values  of  the  resonance  frequencies  from 
those  of  equivalent  electron  systems4.  The  concomitant  increase  in  the 
particle  size  (~  12  A°)  does  not  however  cause  a  large  viscous  broadening  of 
the  plasma  resonance  peaks  since  liquid  helium  below  IK  is  almost  perfectly 
superfluid.  The  observed  width  of  the  resonance  peaks  is  believed  to  arise 
from  interactions  of  the  plasma  waves  with  capillary  waves  on  the  liquid 
surface.  At  typical  experimental  temperatures  of  1  1U  mK,  even  these 
interactions  are  small  and  one  can  easily  obtain  sharp  peaks  with 
Q  -  102  -  103. 

Experiments  have  also  been  done 3  to  study  the  effect  of  a  constanc 
magnetic  field  B0z  perpendicular  to  the  plane  of  the  ions  on  the  wave  proper¬ 
ties  of  the  system.  The  cold  fluid  theory  is  extended  here  to  describe  the 


effect  of  the  magnetic  field  by  including  the  Lorentz  force  term  in  the  equa¬ 
tion  of  motion  for  a  fluid  element.  It  is  interesting  to  note  that  in  the 
presence  of  the  magnetic  field,  the  plasma  equilibrium  can  be  either  static 
or  dynamic.  For  example,  in  the  unmagnetized  cold  plasma  equilibrium 
described  earlier,  each  fluid  element  is  at  rest  since  the  total  electric 

A 

field  acting  on  it  is  zero.  Hence,  imposing  a  magnetic  field  BoZ  does  not 
change  the  equilibrium  density  profile.  On  the  other  hand,  one  can  also  have 
dynamical  equilibria  with  a  magnetic  field.  For  instance,  if  in  the  static 
magnetized  equilibrium,  the  guard  ring  potential  is  changed,  then  each  fluid 
element  would  experience  a  radial  electric  field.  The  plasma  which  would  have 

A 

adjusted  its  radius  in  the  absence  of  BqZ  to  cancel  this  electric  field  now 
executes  an  £  x  ft  rotation  around  the  axis  of  the  cell.  Since  the  radial 
electric  field  is  not  proportional  to  the  distance  from  the  axis,  the  E  *  B 
rotation  velocity  has  a  shear.  Viscosity  effects,  if  included,  can  also  cause 
a  radial  redistribution  of  the  plasma  leading  to  a  more  complicated  equili¬ 
brium  problem.  The  presenc  study  is  limited  to  the  simple  case  of  static 
equilibria. 

The  presence  of  the  magnetic  field  breaks  the  azimuthal  symmetry  of  the 
unmagnetized  modes,  i.e.,  modes  with  azimuthal  mod'  numbers  and  -l 
( l  *  0)  which  are  degenerate  in  the  absence  of  the  magnetic  field  now  have 
different  frequencies.  These  are  the  two-dimensional  analogs  of  modes  found 
in  magnetized  single-component  plasma  columns.6  As  in  the  three-dimensional 
case,  the  lowest  branch  (for  each  value  of  i)  is  an  edge  mode,  whose  existence 
depends  on  the  presence  of  a  free  plasma  edge.  In  the  limit  of  a  large 
magnetic  field,  the  density  perturbation  associated  with  this  mode  is 
localized  at  the  plasma  edge.  The  fluid  motion  is  just  the  E  *  B  motion  and 


so  Che  mode  can  be  regarded  as  a  two-dimensional  analog  of  Che  diocotron  mode. 

The  paper  is  organized  as  follows.  In  Sec.  II,  the  geometry  of  a 
typical  experiment  is  described  and  it  is  shown  that  the  two-dimensional 
approximation  of  the  ion  system  is  justified.  In  Sec.  Ill,  the  horizontal 
equilibrium  of  the  charge  disk  is  obtained  by  balancing  the  self-repulsion  and 
pressure  by  the  confining  electric  fields  produced  by  the  voltages  on  the  cell 
walls.  The  theory  points  to  an  almost  rectangular  density  profile  with  a 
Debye  sheath  near  the  edge.  The  cold  fluid  equations  are  used  in  Sec.  IV  to 
obtain  the  unmagnetized  normal  modes  supported  by  a  plasma  disk  with  a 
rectangular  radial  density  profile.  These  equations  are  extended  to  treat  the 
magnetized  inodes  in  Sec.  V. 


,»4  „y 
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II.  Experimental  Geometry 


A  typical  cell  used  in  the  two-dimensional  ion  plasma  experiments  is 
shown  schematically  in  Fig.  1.  It  consists  of  two  horizontal  electrodes  A  and 
B  a  distance  h  apart,  electrically  insulated  from  a  guard  ring  G  of  radius  R. 
The  cell  is  partially  immersed  in  liquid  4He  (a  screen  grid  at  the  center  of  B 
allowing  the  liquid  into  the  cell)  so  that  the  height  d  of  liquid  helium  is 
roughly  half  the  cell  height.  Helium  ions  are  created  by  a  field  emission  tip 
(positioned  just  below  the  screen  grid  in  B)  which  is  biased  a  large  positive 
voltage  relative  to  B.  Once  the  ions  enter  the  cell  through  the  screen  grid, 
they  are  pushed  up  towards  the  surface  by  the  electric  field  due  to  the 
potential  difference  V^g  between  A  and  B.  However,  the  ions  also  polarize  the 
surrounding  liquid  helium  (dielectric  constant  e  =  1.06)  and  this  polarization 
pushes  the  ions  downwards  with  a  force  which  can  be  simulated  by  image  charges 
(e  -  l)q/(e  +  1)  =<  0. 03q  located  as  far  above  the  liquid  surface  as  the  ions 
lie  below.  The  two  forces  create  a  potential  well  for  the  vertical  motion  of 
the  ions.  The  bottom  of  the  potential  well  is  typically  4  x  10~6  cm  below 
the  liquid  surface  (and  about  20%  further  down  at  the  plasma  edge  where  the 
imposed  vertical  electric  field  is  weaker).  At  typical  temperatures 

3  *  10-2  K,  the  r.m.s.  deviation  of  ions  about  the  bottom  of  the  well  is 
less  than  2  *  10~7  cm.  A  positive  potential  Vq  applied  to  the  guard  ring  G 
provides  the  radial  confinement  of  the  plasma  sheet  which,  typically,  has  a 
radius  ~  1.4  cm  and  contains  ~  6  x  10®  ions  (corresponding  to  an  interparticle 
spacing  of  l  x  10~4  cm).  Thus  the  disk  formed  by  the  ions  is  two-dimensional 
to  an  excellent  approximation. 
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III.  Equilibrium 

The  azlmuthally  symmetric  horizontal  equilibrium  of  the  two-dimensional 
plasma  layer  is  assumed  to  result  from  the  balance  of  the  electrostatic  forces 
on  Che  plasma  and  pressure.  Ignoring  correlation  effects,  the  equilibrium 
charge  density  aQ(r)  thus  depends  on  the  electrostatic  potential  <p0  in  the 
plane  of  the  charges  through  the  Boltzmann  factor  exp(-q  <p0/'[ )  where  q  is  the 
charge  of  a  single  ion  and  T  is  the  temperature  of  the  plasma  (which  is  taken 
to  be  the  temperature  of  the  surrounding  liquid  helium).  The  equilibrium  in 
the  vertical  direction  z  is  not  explicitly  considered  since,  as  indicated  in 
the  previous  section,  the  plasma  layer  is  a  flat  disk  of  virtually  zero 
thickness  located  close  to  the  liquid  surface  for  typical  experimental 
parameters.  Thus,  Poisson's  equation  can  be  written,  to  a  good 
approximation,  in  the  form 

ia  a  a2*o 

7  IF  r  IF  +  =  -4irclao(r)  S(z  -  d) 

*  -4*q<xo(0)  exp  {-  q[$0(r,z=d)  -  <ji0( r=0, z=d )  J /T } 6(z  -  d)  ,  (1) 


! 


% 


where  <(>0  satisfies  the  boundary  conditions  $0  (r  =  R,  z)  =  0, 

♦o  (r,  z  =  0)  =  Vab  -  Vg  and  (r,  z  =  h)  =  -Vq  (where  the  zero  of  the 

potential  has  been  redefined  for  convenience).  It  is  easy  to  verify  that 

CO 

$0  (r,  z)  =  l  [Fn  f(jn  (z)  +  gn  (z)]  JO  UOnt)  ,  (2) 

n=  1 


where 


1  1 

sinh 

k0nd 

sinh 

k0n(h-z), 

for 

z>d 

sinh  k  d  sinh  k  (h-d)  i 
On  (jn  1 

sinh 

R0nz 

si  nti 

R0n(h_d)’ 

for 

z<d 

(3) 


and 


3 
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gn(z) 


2  sinh  k.Qn(z  -  d) 

ko„E  W0 


VG/sinh  k^Ch  ~  d),  for  z>d 


(VAB  '  VG)/Sinh  kUnd’  f°r  Z<d 


,  (*) 


[with  kQn  being  the  solutions  of  Jo(kQnR)  =  0]  is  the  general  solution  of 
Eq.  (1)  for  z  *  d,  which  satisfies  the  boundary  conditions  and  is  continuous 
across  the  layer  z  =  d. 


Integration  of  Eq.  (1)  across  z  =  d  gives  the  jump  condition 


1*2.  I..  "1*2.  L  *  -4*q<*o(r)  , 

3z  'd+  TT  ld~ 


which  together  with  Eqs.  (2)  -  (4)  yields  the  following  nonlinear  equation  for 
the  coefficients  Fn: 

R2  r2  ,  ,  r  k0nsinh  k0nh 

2  i  0nR)  l*n  LSinh  k^d  sinh  k^Ch  ~  d)  ' 


2  ,  VG  VAB  ~  V(j 

R  J  (k_  R)  sinh  k  (h  -  d)  sinh  kA  d 
1  On  On  On 


]}  =  4rrq  /  rdr  J^k^r)  ojc) 


4nq  /  rdr  V^r)  Oq(0)  exp  {-  f  [l  Vo(kOmr)  "  ^  Fm  ^ 
0  mm 


Equation  (6)  can  be  numerically  solved  for  Fn  and  hence  for  <J0(r)  by 
successive  iteration. 

The  numerical  solutions  indicate  that  as  the  temperature  is  decreased, 
the  density  profile  0o(r)  is  almost  constant  up  to  a  radius  (which  is 
determined  by  the  geometry,  the  potentials  applied,  and  the  total  charge) 
where  it  deceases  to  zero  rapidly.  For  d  -  |  h  «  K,  the  scale  length 
for  the  density  fall-off  region  in  the  small  T  limit  is  determined  entirely 
by  h,  the  height  of  the  cell.  Analytic  solutions  obtained  in  the  Appendix 
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for  a  simple  Cartesian  geometry  (in  the  T  -*■  U  limit)  supports  these 
conclusions.  The  analytic  solutions  also  suggest  that  the  T  =  0  density 
profile  near  the  plasma  edge  r  =  a  has  the  form  0o(r  =  a)  ~  (a  -  r)1/2. 

For  non-zero  temperatures,  a  Debye  sheath  [with  Debye  length  defined  as 

*  (Th/4ifq20o(O) ) 1/2]  is  added  on  to  the  plasma  edge.  Figure  2  displays  the 
numerically  determined  equilibrium  profiles  0o(r)  for  X^/h  =  0.1,  0.2  and  for 
typical  experimental  values  VG/oo(0)qh  =  11.6,  VAB/0o(O)qh  =  9.9,  d/h  =  0.467 
and  R/h  =  5;  these  profiles  are  obtained  by  solving  Eq.(6)  on  a  radial  grid 
of  101  points  with  a  basis  set  of  82  Bessel  functions.  The  inner  profile,  for 
which  Xp/h  =0.1,  is  a  good  approximation  to  the  T  =  0  density  profile. 

Typical  experimental  values  of  Xjj/h  are  less  than  3  x  lO-^  and  therefore  the 
plasma  is  described  well  by  the  T  =  0  density  profile. 

In  the  T  =  0  limit,  there  is  a  one-to-one  correspondence  between  the 
values  of  the  central  density  ao(0)  and  the  plasma  radius  a,  for  a  given 
geometry  and  a  given  set  of  confinement  potentials.  This  relation  can  be 
obtained  to  a  good  approximation  by  noting  that  the  T  =  U  profile  is  nearly 
rectangular  and  that  there  is  no  potential  variation  in  the  charge  layer  for 
r  <  a  (i.e.,  no  horizontal  electric  field  for  r  <  a).  For  a  rectangular 
profile  of  radius  a,  one  can  readily  find  the  value  of  0O(U)  needed  to  produce 
a  potential  drop  (due  to  self-electric  field)  between  r  =  0  and  r  =  a  which 
cancels  the  potential  drop  due  to  the  external  confinement  fields.  The 
result  is  a  nearly  exponential  dependence  of  0o(O)  on  a/«  as  displayed  in 
Fig.  3  for  the  values  vAb  =  0.86  VG,  d/R  =  0.093  and  h/K  =  0.2.  The  slope 
of  the  in  Oo(0)  vs  a/R  plot  is  approximately  nR/h.  The  maximum  experimental 
value  of  in [ 0O ( 0)qh/VG J  that  can  be  obtained  experimentally2  is  -2.45 
corresponding  to  the  maximum  possibLe  value  of  the  plasma  radius  a  =  U.895R  = 
1.34  cm.  The  presence  of  a  surface  tension  meniscus  at  the  edge  of  the  liquid 
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helium  pool  is  believed  responsible  for  limiting  the  plasma  disk  radius  to 
to  this  value. 

As  mentioned  in  the  Introduction,  the  static  equilibrium  density  profile 
°o(r)  obtained  in  this  section  is  independent  of  whether  a  magnetic  field  is 
present  or  not.  Thus,  it  can  be  used  to  study  both  unmagnetized  (Sec.  IV)  and 
magnetized  (Sec.  V)  plasma  oscillations  of  the  ion  disk.  Since  the  width  of 
the  density  fall-off  region  is  small  compared  to  the  plasma  radius,  the 
equilibrium  density  profile  can  be  approximated  well  by  a  rectangular  profile. 
This  simplifying  assumption  is  used  in  the  following  sections. 
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V.  WAVES  IN  THE  ABSENCE  OF  A  MAGNETIC  FIELD 

A.  Theory 

Since  the  equilibrium  is  stationary  in  time  and  is  azimuthally  symmetric, 
the  perturbed  density,  velocity,  and  potential  can  be  assumed  to  be  propor¬ 
tional  to  exp(-ioit  +  i£9)  and  thus  the  linearized  continuity  equation, 
equation  of  motion  and  Poisson's  equation  take  the  form 


-icon.  +  —  4—  (ru  v,  J  +  o  —  v .  0  =  0 
1  r  3r  o  lr  or  10 


(7) 


-ioiv 
~  11 


(8) 


r32  +  j_  3_  3 _ iK 

3z2  r  3r  3r  r2 


—4  TTq  a 


5(  z 


-  d)  , 


(9) 


where  CJQ ,  <Jj  and  =  v^r  +  are  functions  of  r,  and  4>j(t,z)  is  the 

mode  potential  and  $e(r,z)  exp( -io)t+i  £9)  is  the  potential  due  to  the  driver 
voltage  applied  to  the  walls.  It  should  be  noted  that  Eqs.  (7)  and  (8)  are 
two-dimensional  equations  defined  only  in  the  plane  (z  =  d)  of  the  charge 
layer,  while  Eq .  (9)  is  a  three-dimensional  equation.  To  complete  the  problem, 
boundary  conditions  on  have  to  be  specified.  These  are  obtained  by  noting 
that  the  walls  of  the  cell  are  conductors  and  that  the  impedance  between 
different  sections  of  the  wall  at  typical  mode  frequencies  is  small;  thus 
can  be  assumed  to  vanish  at  the  boundaries  z  =  0,  z  =  h,  r  =  R  of  the  cell. 


For  z  *  d,  Eq .  (9)  reduces  to  Laplace's  equation  for  and  the 
general  solution  satisfying  the  boundary  conditions  and  being  continuous  at 
z  *  d  is 


:w*ivi*Aim***n**i***i**p»jruwtf  mj  w*u»uwv^\rr\<Tvrv?njirvjwvw\/*virvrvw;iP.hFJWJ»« 
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*l(r,z)  ~Jn  A"  f*n(z)  J*(k*nr) 


(1U) 


where  k  are  the  roots  of  J  (k  R)  =  0  and 
in  i  tn 


f  (z)  »  _ i _ 

in'  sinh  k.  d  sinh  k.  (h 


>-  fk  U  0  1(111  cv  -  \ 

Xn  Xn 


sinh  k,  d  sinh  k  „  (h 
in  in 


sinh  k.  z  sinh  k.  (h 
in  in 


z) ,  z  >  d 
d) ,  z  <  d 


(ID 


is  defined  such  that  f  (z~d)  =  1.  The  functions  Oi(r)  and  <f»e(r,  z 

in 


d)  are 


also  expanded  in  terms  of  J  (k  r): 

i  in 


a  (r )  =  l  B  J  (k  r)  , 
1  n* 1  n  i  in 


(12) 


$  (r,  z  =  d)  =  l  C  J  (k  r)  . 
e  n= l  n  i  in 


(13) 


Integrating  Eq.  (9)  across  z  =  d  gives  the  jump  condition 


9$. 


3z 


3*. 


d+ 


3z 


d- 


-4wqo1(r ) 


(U) 


Combining  Eqs.  (7)  -  (14)  yields  the  relation 
A 


i  =  -  l  (D . “1 )  (D-)  C 

n  L  1  nm  2  mp  c 

m,p 


(15) 


between  the  Four ier-Besse 1  components  A^  0f  ^  (r>z=d)  and  the  components 
Cn  of  the  known  driver  potential  $e(r,z=d),  with  the  matrices  Ui  and  D<> 
defined  by 


(U,  ) _ =  (D.) 


?  k  „  R  sinh  k„  h 

U)^  R  in  in 


-  v-o/  „  ,  ,  -  J  (k  R)  6  (lb) 

1  mn  2  mn  2oip 2  h  sinh  kSnd  sinh  k>ri(h  -  d)  i+1  in  mn 


‘ink¬ 


in' 
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In  these  definitions,  6mn  is  the  Kronecker  delta  function,  primes  repre¬ 
sent  derivatives  with  respect  to  the  arguments,  and  oD  are  defined  as 

„  4irq2o  (0) 


.  <3  (r) 

°o(r)  *  ffTJTT  *  (19) 

o 

Equations  (15)  and  (16)  imply  that  the  plasma  response  has  resonances  at 

frequencies  for  which  Det  (Dj)  =  0.  The  resonance  frequencies,  when  normalized 

in  terms  of  are  completely  specified  by  the  values  of  h/R,  d/R  and  aQ  (r/R). 
The  numerical  solutions  for  the  lowest  few  resonance  frequencies  when  oQ(r)  is 

assumed  to  be  a  step  function  of  radius  a  are  displayed  in  Fig.  4  for  typical 

experimental  values,  h/R  =  0.2,  d/R  =  0.093,  a/R  =  0.895.  Each  plasma 

resonance  is  characterized  by  an  azimuthal  mode  number  i  and  a  radial  mode 

number  N.  The  (numerically  smoothed)  functions  <frl(r,z=d),  o  (r),  v  (r)  and 

1  Ir 

v  (r)  are  displayed  in  Fig.  5  for  a  frequency  close  to  the  l  =  1,  N  =  2 
1  0 

resonance  chosen  as  a  typical  example.  Since  from  Eqs.  (7)  and  (8),  oi 
contains  a  term  proportional  to  3o0/ 3r *3$i/ 3r  =  -ao(0)6(r  -  a)3$l/3r  and  since 
the  radial  electric  field  does  not  vanish  at  r  =  a,  has  a  sharp  spike  at  the 
plasma  edge.  It  is,  however,  noticed  in  the  numerical  solutions  (as  in  Fig.  5) 
that  $^(r<a,z=d)  has  the  form  of  a  single  Bessel  function  with  a  rather  small 
derivative  at  r  =  a  (as  long  as  a/R  is  not  very  close  to  unity).  In  fact,  the 
assumption  =  ANJJl(K^Nr)f  ^N(z )  [where  is  a  solution  of  J^(Ka)  =  0  ai  J 

A  » 

fgN  is  defined  by  Eq .  (11)  with  taking  the  place  of  k^]  can  be  used 

in  Eqs.  (7)  -  (9)  to  obtain  good  approximations  to  the  resonance  frequencies 
of  a  step-function  equilibrium  density  profile: 
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K^Na  sinh  K^d  sinh  K^N  (h  -  d) 
sinh  K'  0l,  h 


(20) 


Equation  (20)  predicts  that  when  plotted  versus  Kj^a,  the  resonance 
frequencies  for  all  t  fall  on  the  same  curve.  The  solid  line  in  Fig.  6 

represents  this  curve  for  d/a  =  0.104,  h/a  =  0.223.  The  exact  numerical 
solutions  (d/a  =  0.106,  h/a  **  0.223,  a/R  =  0.895)  using  the  expansion  (10) 
for  lie  very  close  to  this  curve.  The  assumption  d$/dr  =  0  at  the  con¬ 
ducting  wall  (r  =  R)  has  been  used  previously  ^  to  obtain  an  equation 
similar  to  Eq.  (20)  which  agrees  well  with  the  experimental  values  of 
resonance  frequencies.  Figure  6  shows  why  this  assumption  works  well  (even 
though  it  is  unjustified)  when  a  is  nearly  equal  to  R.  The  limits  of 
applicability  of  Eq.(20)  are  discussed  later. 

The  simple  approximate  structure  of  suggested  by  the  numerical 
results  can  be  understood  by  considering  Eqs.  (7)  -(9),  which  can  be 
combined,  ignoring  4»e ,  as 

u.2V2(4>  e1*9)  =  V  •  a  (r)  V  (♦  eiW)]6(z  -  d)  .  (21) 

1  X  m  o  1  1 


Multiplying  both  sides  by  <)>jexp(-i  £9)  and  integrating  over  the  cell  volume 
yields  the  variational  (Rayleigh-Ritz )  principle  value  of  u)2  , 
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The  minimum  value  of  the  right-hand  side  gives  the  lowest  eigenvalue  u>2  and 

the  $1  that  makes  it  the  minimum  is  the  corresponding  eigenfunction.  The 

second  lowest  eigenvalue  is  the  minimum  produced  by  a  wave  function  which  is 

orthogonal  to  the  first  and  so  on.  It  is  of  interest  to  note  that  Eq.  (22) 

falong  with  Eq.  (8)]  implies  that  the  total  energy  in  the  wave  electric  field 

is  equal  to  the  total  kinetic  energy  of  the  plasma  particles. 

To  proceed  further,  the  denominator  in  Eq.  (22)  is  converted,  on 

integration  by  parts,  into  the  area  integral  4itq  ^  rdr  $1  (r,z  *  d). 

o  ’ 

From  Eqs.  (14),  (11),  (10),  and  (12),  it  follows  that  if  k  h  «  1,  then 

in 


D  2  »  » 

n  4irq  d(h  -  d)  n  ’ 

thus  if  h/R  is  small  and  if  higher  order  Fourier-Bessel  components  can  be 
ignored,  then 


oL(r) 


1  h 
4irq  d(h  -  d) 


z=d  ) . 


(23) 


Thus  Eq.  (22)  takes  the  form 


/ 

o 


a)2  =  ui  2  d(h  -  d) 

P  a 


7l^l^r,Z=d^ei*9^ 


rdr 


$^(r ,z=d)e 


iZ9 


rdr 


(24) 


The  minima  of  the  ratio  of  the  integrals  can  be  verified  to  be  [which 

occur  for  in  a8reement  with  Eq.  (20)  in  the  limit  of  small 

h/R. 
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The  discussion  above  suggests  that  Eq.  (20)  might  fail  when  the  value  of 
h/R  is  not  small  compared  to  unity.  Figure  7  displays  how  the  deviation 
between  the  prediction  of  Eq.  (20)  (solid  line]  for  the  Z  =  0,  N  =  2  mode 
chosen  as  an  example  and  the  exact  numerical  solutions  of  Eqs.  (15)  -  (17) 
[circles]  increases  with  increasing  values  of  h/R.  Figure  8  presents  the 
comparison  between  the  predictions  of  Eq.  (20)  [solid  line]  for  the  same  mode 
and  the  exact  numerical  solutions  [circles]  as  a  function  of  a/R  for  h/R  «  0.2. 
The  percentage  difference  is  relatively  constant  until  a/R  *  0.98,  beyond 
which  it  rapidly  increases.  The  limiting  value  of  the  exact  solution  as  the 
plasma  edge  approaches  the  walls  can  be  obtained  from  Eqs.  (15)  -  (17).  For 
a/R  ♦  1,  the  matrix  D2  becomes  diagonal, 

2  2 

(Do)  «  Jl+l(kJlmR)  ^mn  .  (25) 

‘'mn  2 


and  thus  the  resonance  frequencies  (which  are  the  zeros  of  D^)  have  the  values 


,  „  h  k*N  R  Slnh  k*N  d  Slnh  k*N(h'd) 

4N  P  R  sinh  k.^  h 


(2b) 


corresponding  to  wave  functions  6  =  B  J  (k  r)f  (z)  which  vanish  at 

1  N  l  iN  *N 

r  »  R.  It  is  found  numerically  that  the  transition  from  the  modes  of  Eq. 

(20)  (whose  radial  derivative  at  the  plasma  edge  is  small)  to  the  modes  of  Eq , 
(26)  (whose  value  at  the  plasma  edge  vanishes)  occurs  for  values  of  a/R  very 
close  to  unity  (a/R  >  0«98  for  typical  values  of  h/R).  The  frequencies  of  Eq . 


(26)  are  roughly  (k  /K^)^  higher  than  the  values  of  Eq .  (20).  Since  in 


experiments,  the  plasma  radius  a  cannot  approach  the  cell  radius  R  to 
within  M.5  mm  due  to  the  liquid  helium  meniscus,  modes  described  by  Eq .  (2b) 
cannot  be  observed  in  cells  with  R  <  8  cms. 


L 
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B.  Comparison  with  Experiments 

Waves  are  excited  in  the  experiments  of  Ref.  1  by  imposing  an  oscillating 
potential  on  the  guard  ring  and  are  detected  by  measuring  the  current  from 
a  circular  button  (of  radius  rc)  at  the  center  of  the  top  electrode.  If 
the  apparatus  is  exactly  cylindrically  symmetric,  only  l  *  U  modes  should  be 
detected.  The  values  predicted  by  Eq.  (20)  for  the  first  four  radial  modes 
having  no  angular  dependence  agree  with  the  experimental  values  co  within  a 
few  percent.  Other  resonances  have  also  been  observed2,  which  can  be 
identified  with  the  predictions  of  Eq.  (20)  for  various  low  values  of  l  {*  0) 
and  N;  these  probably  arise  from  an  undetermined  azimuthal  asymmetry  in  the 
system. 

For  the  1*0  modes,  the  absolute  value  of  the  current  flowing  out  of  the 
center  button  when  the  driver  frequency  matches  a  resonant  frequency  [given 
approximately  by  Eq.  (20)],  can  be  obtained  by  including  the  collisional 
damping  term  v^li  on  the  left-hand  side  of  Eq .  (8).  This  has  the  effect  of 
replacing  w2  in  Eq.  (16)  by  w(w  +  iv).  If  v  is  small  compared  to  the 
frequencies  w  of  interest,  this  is  equivalent  to  replacing  w  by  w  +  iv/2  in 
Eq.  (16).  For  w  near  won  [which  are  simple  zeros  of  Di(w)  as  follows  from  the 
self-adjointness  of  Eq.  (21)],  one  can  write 

°i“i  (a> +  iv/2)  2  —  'o^r;  Tw 2  url  (uluN  +  v)  • 

where  Dj-1  on  the  right  hand  side  can  be  conveniently  evaluated  at  a  real 
value  Won  +  v  of  Che  frequency.  From  Eq .  (15)  a  similar  relation  holds 
for  the  coefficients  An  evaluated  at  w  +  iv/2.  Therefore,  from  Eq.  (10) 


the  wave  potential  has  the  form 
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#l(r’  z)  "Jj  I"lv/2  foN  (2)  J°  (k0Nr) 


when  u>  Is  close  to  o»qN  and  v  is  snail  compared  to  cd.  The  real  coefficients 
+  are  obtained  by  numerically  solving  Eq.  (15)  for  id  =  idqn  +  v 


^n  *  -2<fry  [sinh  k  d  +  sinh  k  (h-d)j  [k  RJ[  (k  R)  sinh  k  h]“  . 

On  On  On  On  On 

The  charge  density  induced  on  the  center  button  by  the  plasma  wave  potential 
is  q^ind  3  (4w)_1 3$i/3z | 2»h  dnd  the  magnitude  of  the  current  through  the 
center  button  is  given  by 


I(u>  a  (DUN)  *  <d[ 2w  /  r  dr  |otnd|] 
0 


v  An  (<DQN  +  v) 


2  rc)  lmi  [(„  -  <dqN) 2  +  v2/4J 


Jl  (k0Nrc) 
sinh  kuN  (h-d) 


For  cd  not  too  close  to  <*>qn,  one  can  ignore  iv/2  and  replace  the  terms  in  the 
braces  by  A^fo)).  The  plot  of  I(<d)  vs.  <d  consists  of  a  series  of  Lorentzian 
resonance  peaks  centered  at  cdqN  and  of  width  v.  This  is  displayed  in  Fig.  9 
for  a  typical  value  of  the  collision  frequency  v  =  4  x  10-3  ^  and  for  a  range 
of  id  covering  the  first  four  l  *  0  resonances. 


. .  • t.f 
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V.  WAVES  IN  THE  PRESENCE  OF  A  MAGNETIC  FIELD  ft  -  Bqz 


Imposing  a  uniform  magnetic  field  ft  a  Bqz  perpendicular  to  the  plane  of 
the  plasma  disk  introduces  the  Lorentz  force  term  Q  x  2  (where  £1  ■  qB0/mc 
is  the  cyclotron  frequency)  to  the  right-hand  side  of  the  linearized  equation 
of  motion  (8).  Thus  the  perturbed  velocity  components  are  given  by 


'  i  _  “  — s - — y  (—iw  ~ —  +  — — — )  [♦,+♦]  , 

ir  u)2  -  ft*  m  9r  r  1  e  z*d 


V10  "  ~  2  ~~  ^  (~  _  R  f"5  +  ♦  1  „ 

nr  *  U2  a  r  3r  1  e  z-d 


(29a) 


(29b) 


Using  the  expansions  (10)  -  (13)  for  <fri,  a\,  and  $e,  one  again  obtains  Eq. 
(15),  but  with  the  matrices  Dj  and  D2  now  given  by 


2  r>2  a  k.  R  sinh  k,  h  ~ 

(D  )  -  (D  )  -  — -~-Q-  -  - — - ^ -  J2  (k  R)  5  (30) 

1  mn  2  mn  2u»p2  h  sinh  k^d  sinh  kjn(h  -  d)  i+1  in  mn  , 


<D2)mn  *  l  rdr  5o(r)  [kimkin  Ji(kimr)  JiUinr)  +  jz  Vktn»r)  V^01 


i£2  k  d0o(r) 
+  “  l  rdr  ~dr~ 


Ji(kimr)  Ji(kinr) 


If  i  *  0  or  if  Oo  is  constant  for  0  _<  r  R,  then  the  last  term  of  Eq .  (31) 
vanishes;  then  Eqs.  (30)  and  (31)  can  be  obtained  from  Eqs.  (16)  and  (17)  by 
replacing  w2  by  u>2  -  ft2.  Thus,  the  square  of  the  resonance  frequencies 
in  the  magnetized  case  are  shifted  up  by  an  amount  n2  from  the  unmagnetized 
values  with  the  eigenfunctions  remaining  unchanged.  These  are  the 


two-dimensional  analogs  of  the  usual  upper  hybrid  resonance. 


w 


Without  the  last  term  in  Eq.  (31),  modes  corresponding  to  +  A  and  -1 l  are 
degenerate.  For  t  *  0  and  dSQ/dr  *  0  (for  some  r<R),  the  last  term  in  Eq. 

(31)  breaks  the  degeneracy  and  causes  each  unmagnetized  mode  to  split  into 

two  upper  hybrid  modes  which  propagate  in  opposite  directions  with  azimuthal 

phase  velocities  tv/ A.  In  addition,  two  modes  appear  near  a >  “  52,  corresponding 

to  +A  and  -A.  The  effect  of  the  magnetic  field  for  | A|  ■  1  is  shown  in  Fig. 

10  for  a/R  *  0.895,  d/R  =  0.093  and  h/R  *  0.2,  where  a  is  the  width  of  the 

A 

rectangular  density  profile  0o(r).  The  modes  occur  in  pairs  designated 
+  or  -  depending  on  whether  they  they  propagate  circularly  right  handed  or  left 
handed  with  respect  to  the  magnetic  field.  The  two  new  modes  near  <v  =  52  are 
also  shown.  It  should  be  noted  that  while  these  modes  occur  very  close  to  the 
cyclotron  frequency  52,  there  is  no  resonance  when  tv  =  52  exactly.  It  is  of 
interest  to  note  that  when  tv  =  52  exactly,  Eqs.  (30)  and  (15)  yield 
and  An  »  -Cn.  Thus  the  density  perturbation  produced  in  the  plasma  is  such 
chat  the  total  wave  potential  +  <fre  is  zero  in  the  plane  of  the  charges. 

This  is  true  for  all  values  of  A  including  A  =  0. 

The  forms  of  <t>i(r,  z  =  d)  and  0j(r)  for  the  magnetized  modes  are 
qualitatively  similar  to  the  corresponding  magnetized  modes.  The  wave 
potential  $i,(r,z  *  d)  for  the  positive  A  mode  near  <v  =  52  increases  linearly 
with  r  till  the  plasma  edge  and  then  decreases  to  zero  at  r  =  R;  the  perturbed 
density  0j  increases  monotonically  to  a  peak  at  r  =  a,  as  displayed  in  Fig.  11 
for  the  case  A  *  +  1.  The  curves  are  obtained  by  numerically  solving  Eq .  (15) 
for  h/R  =  0.2,  d/R  =  0.093,  a/R  =  0.895,  52/ (Vp  =  n.5  and  u>/ u>p  =  0.504.  A  basis 
set  of  82  Bessel  functions  is  used  on  a  101-point  radial  grid;  the 
oscillations  in  for  r  >  a  and  the  finite  height  of  0[  at  r  =  a  are  con¬ 

sequences  of  these  numerical  limitations.  For  the  negative  A  mode  near  tv  =  51, 
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Che  wave  potential  ^[(r,  z  =  d)  is  mainly  confined  to  the  annulus  a  <  r  <  R 
and  the  perturbed  density  0^  is  localized  at  the  plasma  edge  r  =  a. 

The  behavior  of  the  magnetized  modes  is  similar  to  the  unmagnetized  modes 
as  a/R  approaches  unity.  In  this  limit,  the  last  term  in  Eq.  (31)  vanishes, 
the  matrix  D2  becomes  diagonal,  and  the  resonance  frequencies  take  the  values 


01 
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°2  -  V  if 


k^R  sinh  k^d  sinh  k^N(h-d) 
sinh  k^N  h 


(32) 


As  in  the  unmagnetized  case,  the  frequencies  of  Eq.  (32)  correspond  to 

eigenmodes  (r,  z  *  d)  ~  J  (k  r).  It  is  worth  noting  that  without  a  free 
1  l  4N 

plasma  edge,  the  +4  and  -t  modes  are  degenerate  even  in  the  presence  of  a 
magnetic  field. 

As  a/R  increases  towards  1,  the  frequencies  of  the  i  l  modes  (for  any 

given  value  of  £1)  shift  to  the  values  oj  given  by  Eq.  (32)  accompanied  by 

JIN 

a  corresponding  change  in  the  form  of  $i(r,z  *  d);  the  exceptions  are  the  -l 
mode  near  ui  =  Q  which  approaches  u>  =  Q,  and  the  lowest  +  l  mode  which 
approaches  ui  =  0.  As  in  the  unmagnetized  case,  the  change  is  rather 
abrupt  near  a/R  =  0.99.  As  shown  in  Fig.  10  (for  | l\  =  1),  the  bottom  two  +1 
modes  cross  each  other  and  it  becomes  clear  upon  letting  a/R  approach  unity 
that  there  is  a  mode  transition  at  the  crossing.  This  is  illustrated  in  Fig. 
12  for  the  case  l  =1.  The  two  upper  lines  go  over  to  u)j  j  given  by  Eq.  (32) 
in  the  limit  a/R  •+•  1  and  must  be  regarded  as  the  same  mode.  Similarly  the 
lower  two  lines  approach  u>  =  0  in  the  limit  a/R  =  1  and  hence  belong  to  the 
same  branch.  For  the  lowest  branch,  for  all  l  *  0,  the  ratio  u>/ £2  approaches 

zero  for  large  values  of  ft  and  in  this  limit,  Eq .  (29)  yields 

”  “  (c/B)  Vj_  ($j  -t-  $e)z=(j  x  z  indicating  incompressible  motion 


(7  ,  v  *0)  along  the  plane  of  the  charges.  Also  from  the  continuity  equa- 
X  ^11 

tion,  *  (ioj)-i  v^r  dcQ/dr  =  -(ioi)~l  v^r  ao(0)  6(r-a);  thus  the  perturbed 
charge  density  is  localized  at  the  plasma  edge  r  =  a.  The  mode  in  the  limit 
of  large  is  purely  an  edge  mode  and  can  be  interpreted  as  a  two-dimensional 
analog  of  the  usual  diocotron  mode.  It  should  however  be  noted  that  the 
two-dimensional  equilibrium  considered  here  does  not  execute  an  E  x  B 
rotation  and  hence  phenomena  such  as  negative  energy  waves  and  resonance 
layers  are  absent  in  the  present  system. 


VI.  CONCLUSIONS 


A  theory  of  the  equilibrium  and  wave  properties  of  a  two-dimensional  ion 
plasma  confined  just  below  the  surface  of  liquid  helium  has  been  developed. 

In  equilibrium,  it  is  found  that  for  the  typical  confinement  geometry  used  in 
experiments,  the  plasma  density  profile  is  almost  rectangular  in  the  cold 
limit;  for  T  *■  0,  a  Debye  sheath  is  introduced  at  the  edge.  A  theory  of 
linear  waves  that  considers  the  plasma  as  a  cold  two-dimensional  fluid  has 
been  constructed.  Imposing  the  condition  that  the  linearized  wave  potential 
vanish  at  the  conducting  walls  of  the  confining  cell  leads  to  an  infinite  set 

of  radial  modes  for  each  value  of  the  azimuthal  mode  number  l.  Through  a 

variational  approach,  it  is  shown  that  an  excellent  analytical  approximation 
for  these  frequencies  can  be  obtained  (for  h/R  << 1  and  a/R  not  very  close  to 

unity)  and  that  the  wave  potential  in  the  plane  of  the  charges  for  any 

resonant  frequency  has  the  form  of  a  Bessel  function  with  a  rather  small 
derivative  at  the  plasma  edge.  Imposition  of  a  uniform,  constant  magnetic 
field  perpendicular  to  the  plane  of  the  charges  dot  '  not  change  the 
equilibrium  of  the  system  in  the  T  =  0  limit.  However,  the  frequencies  of  the 
1  =  0  normal  modes  are  shifted  upwards  and  for  |  £ |  *  l),  each  normal  mode  is 

split  into  two  modes  (corresponding  to  +  £  3nd  -£)  which  propagate  azimuthally 
in  opposite  directions.  Two  new  modes  (  +  £  and  -l)  appear  for  frequencies  very 
close  to  the  cyclotron  frequency.  The  two  lowest  +£  modes  cross  each  other 
and  it  is  shown  that  there  is  a  mode  transition  at  that  point.  The  lower 
branch  is  an  edge  mode  and  is  identified  as  the  two-dimensional  analog  ot  the 
diocotron  mode.  The  other  branches  are  the  analogs  of  the  familiar  upper 


hybrid  mode 
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APPENDIX 


EQUILIBRIUM  PROFILE  IN  CARTESIAN  GEOMETRY 

A  simplified  model  geometry  in  which  the  equilibrium  density  profile  of 
the  two-dimensional  plasma  can  be  analytically  obtained  in  the  zero 
temperature  limic  is  described  here.  The  model,  shown  in  Fig.  13,  is  a 
Cartesian  model,  one  in  which  the  plasma  is  a  ribbon  infinitely  long  in  the 
direction  perpendicular  to  the  plane  of  the  figure.  The  notation  is  slightly 
different  from  that  of  Fig.  1.  The  vertical  coordinate  is  now  denoted  by  y 
(with  z  used  to  donote  the  complex  variable  x  +  iy).  Also,  for  convenience 
the  plasma  layer  is  assumed  to  be  midway  between  the  top  and  the  bottom  of  the 
confining  box.  As  in  Sec.  II,  the  vertical  equilibrium  is  not  explicitly 
considered  and  consequently  the  vertical  holding  field  is  ignored.  Thus,  the 
model  consists  of  a  ribbon  of  charge  constrained  (by  forces  not  included  in 
the  analysis)  to  the  midplane  of  the  containing  cell.  The  self-repulsion 

of  the  plasma  particles  balances  the  confining  force  due  to  the  potential  VQ 

applied  to  the  side  walls  and  thus  determines  the  horizontal  equilibrium 
density  profile  <Jo(x),  Since  the  temperature  is  assumed  to  be  zero,  the 
plasma  pressure  plays  no  role  in  the  equilibrium. 

Because  the  horizontal  electric  field  on  each  plasma  parcicle  is  zero, 
4>(x,y  =  U)  =  Vj  =  constant  for  |x|  <  a.  From  symmetry,  the  lines  AB  and  UD  in 

Fig.  13  are  lines  of  force.  One  thus  has  the  problem  of  finding  the  charge 
distribution  on  UA  when  it  is  at  potential  the  side  BC  is  at  0,  the 

side  CD  at  VQ  with  AB  and  OD  being  lines  of  force.  Conformal  mapping  is  used 
to  solve  tor  the  potential  distribution  and  hence  the  charge  density  on  UA. 

The  potential  is  not  known  a  priori,  but  it  is  found  that  for  only  one 

particular  value  of  Vj  is  the  calculated  charge  density  finite  everywhere  in 

OA  and  has  the  correct  sign;  this  special  value  of  potential  must  be  that 
corresponding  to  the  equilibrium  in  the  cold  limit. 
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The  interior  of  the  rectangle  OBCD  can  be  mapped  into  the  upper  half  of 
the  complez  5 -plane  by  the  elliptic  function  transformation 


C  •  C  +  in  ■  ra  sn2(K(m)(l  +  > 

d 

where  z  *  x  +  iy.  Defining  the  parameter  m  by  the  equation 

d  K(m) 

T  *  K'(m)  * 


(Al) 


(A2) 


it  is  seen  that  the  points  D,  0,  A,  B,  C  in  the  z-plane  are  mapped  into 

CD  =  0,  Co  “  m,  Ca»  Cb  *  1.  Cc  ”  ®  of  the  5-plane,  with  the  interior  of  the 
rectangle  DOBC  mapped  into  the  upper  half  of  the  5-plane  (Fig.  14a).  To  keep 
the  solution  reasonably  simple,  it  is  also  assumed  that  the  plasma  half-width 
a  is  such  that  Ca  “  1  -  m;  from  Eq.  (Al),  this  implies 

?A  =  m  nd2[K'(m)  1  -  m]  *  1  -  m.  (A3) 


The  problem  thus  reduces  to  finding  a  function  P  =*  +  i$  analytic  in  the 

upper  half  of  the  5-plane  with  $  satisfying  the  conditions  indicated  in  Fig. 
14a,  along  the  real  axis.  The  solution  to  the  problem  is  the  superposition 
of  the  complex  potentials  satisfying  the  boundary  conditions  of  Figs.  14b  and 
14c.  The  respective  solutions  are 
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The  surface  charge  density  of  the  plasma  sheet  is  given  by 


1  d$  =  -1 

2irq  dy  y*0,  |x|<a  2xq 


y=0,  |xj<a 


K’(m) 


1 


(1  -  m  -  5)172  l2K’([l  -  2m]2) 

The  surface  charge  density  has  a  singularity  at  £  =  1  -  m  [or  from  Kq 
x  “  a]  unless 
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in  which  case  the  charge  density  is 
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The  dependence  of  on  x  is  displayed  in  Fig.  15  for  the  value 

d/L  =  0.625  (or  m  =  .1).  To  demonstrate  that  the  warm  plasma  density  profile 
approaches  this  shape  in  the  limit  T  -*■  0,  the  numerical  solutions  of  the 
Cartesian  analog  of  Eqs.  (1)  -  (6)  are  also  shown  in  Fig.  15  for  the  same 
geometry,  the  same  <Jo(0)  and  Xp/d  =  0.05  and  0.1  (where 
XD  -{(2dT/4xq2Oo(0)}1/2]. 


For  small  values  of  m  [or  from  Eq.  (A2),  small  values  of  d/L] ,  Eq.  (A8) 
can  be  expanded  near  x  =  a  to  find  the  behavior  of  the  density  profile  near 


the  edge.  This  yields 


,»i/.  ^+(f 


(A9) 


Thus,  very  near  the  plasma  edge  x  =  a,  oQ  (x)  ~  (a  -  x)1/2.  The  width  of  the 
plasma  edge  is  approximately  2d/ir. 


The  (a  -  x)*/2  dependence  for  oQ(x  3  a)  can  also  be  obtained  from  a 
simpler  argument.  The  solution  of  the  two-dimensional  Laplace's  equation 
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The  charge  density  for  x  <  a  is 


o  (x  £  a)  =  7 
o  p  39 


0*ir  2  X1 


4  Q  (a  -  x)1/2  -  4  Q.  (a  -  x)3/2  +  .  . 


2  2 


(A13) 


in  agreement  with  Eq.  (A9). 


In  the  cylindrical  case,  Laplace's  equation  for  the  equilibrium  potential 
is  (with  z  *  vertical  coordinate  in  the  notation  of  Sec.  Ill) 

+  +  (au) 

3r2  r  3r  3Z2 

which  can  be  locally  (for  r  *  a)  transformed  into  the  form  (A10)  by  setting 
r  *  a  exp(eu),  z  =  eav  in  the  limit  e  >  0.  Following  the  same  procedure, 
one  can  show  that  even  in  the  cylindrical  case,  the  equilibrium  density 
profile  <J0(r  a  a)  *  (a  -  r)1/2. 


*  j  j  y 


i 
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FIGURE  CAPTIONS 


Typical  experimental  geometry. 

Equilibrium  density  profiles  oQ(r)  for  Xy/h  =  0.1,  0.2  and  for  typical 
experimental  values  Vy/oo(0)qh  =  11.6,  VAB/ o0( 0)qh  =  9.9, 
d/h  *  0.467  and  R/h  =  5. 

Central  density  <Jo(0)  required  to  produce  a  plasma  disk  of  radius  a 
for  T  -  0;  VAB  =  0.86Vg,  d/R  =  0.093  and  h/R  =  0.2. 

Unmagnetized  normal  mode  frequencies  for  d/R  =  0.093,  h/R  *  0.2  and  a/k  = 
0.895.  The  modes  are  characterized  by  azimuthal  mode  number  l  and 
radial  mode  number  N. 

Radial  dependence  of  wave  potential  (r,  z=d),  density  a^(r), 

radial  velocity  v  (r)  and  azimuthal  velocity  v  (r)  for  the  4=1, 
lr  19 

N  =  2  mode;  the  geometry  parameters  are  d/R  =  0.093,  h/R  =  0.2  and 

a/R  =  0.895.  If  <b  is  real,  a  and  v  are  also  real  and  v 

l  l  19  lr 

is  purely  imaginary  [see  Eqs.  (7)— (9) ] . 

Comparison  of  the  exact  values  of  the  resonance  frequencies 

and  the  approximate  values  given  by  the  analytical  expression  (20) 

[solid  line]  for  d/a  =  0.104,  h/a  =  0.223  and  a/R  =  0.895. 

Comparison,  as  a  function  of  h/R,  of  exact  values  of  resonance  frequencies 

(circles)  and  the  approximate  values  given  by  the  analytic  expression 
(20)  [solid  line]  for  the  4=0,  N  =  2  mode  and  for  a/R  =  0.895, 
d/h  =  0.465. 

Comparison,  as  a  function  of  a/R,  of  exact  values  of  resonance  frequencies 

(circles)  and  the  values  given  by  Eq .  (20)  [solid  line]  for  the  l  =  0, 

N  =  2  mode  and  for  d/R  =  0.093,  h/R  =  0.2. 
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9.  Current  I  (in  amperes)  through  the  center  button  in  the  top  electrode 
as  a  function  of  the  external  frequency  u»  [Eq.  28].  The  collision 
frequency  v  is  4  x  10~3  uip  and  the  driver  voltage  =  lOmV. 

10.  Plasma  resonances  vs.  ft  =  qB/mc  for  £  *  ±1  and  for  d/r  =  0.093,  h/R  = 
0.2,  a/R  *  0.895.  The  branches  designated  +  (or  -  )  propagate 
circularly  right-handed  (or  left-handed)  with  respect  to  the  magnetic 
field. 

11.  Radial  dependence  of  the  potential  4>i(r,z=d)  and  density  <xi(r) 
for  the  l  *  +1  mode  near  a>  =  ft  for  ft/uip  =  0.5,  d/R  =  0.093, 
h/R  *  0.2,  a/R  =  0.895. 

12.  Enlargement  of  the  region  where  the  l  =  +1  modes  cross  in  Fig.  10. 

The  positions  of  the  branches  as  a/R  assumes  the  values  0.895,  0.985 
and  1.0  are  shown. 

13.  Cartesian  model  geometry  for  the  confined  two-dimensional  plasma  layer. 

14.  a)  Complex  C -plane  showing  the  mapping  (Al).  The  shaded  area  (rectangle 
OBCD)  of  Fig.  (13)  is  mapped  into  the  upper  half  of  the  £-plane.  The 
problem  can  be  regarded  as  a  superposition  of  two  potential  problems 
illustrated  in  (b)  and  (c). 

15.  Density  profile  0o(x)  given  by  the  analytic  expression  (A8)  lor 

d/L  =  0.625  or  m  =  0. 1  (thick  line).  The  other  curves  are  numerical 

solutions  of  the  Cartesian  analogs  of  Eqs.  (1)  -(b)  for  the  same 

geometry,  the  same  0Q( o)  and  for  A^/d  =  0.U5  and  0.1;  the  Debye 

1/2 

length  Au  is  defined  as  Aq  =  [ 2dT/4xq 2oQ( 0 ) ] 

16.  Geometry  in  the  vicinity  of  the  plasma  edge  x  =  a. 
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